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Abstract. 

We present a scheme that produces a strong U(l)-like gauge field on cold atoms 
confined in a two-dimensional square optical lattice. Our proposal relies on two 
£3 ' essential features, a long-lived metastable excited state that exists for alkaline-earth 

or Ytterbium atoms, and an optical superlattice. As in the proposal by Jaksch and 
Zoller [New Journal of Physics 5, 56 (2003)], laser-assisted tunneling between adjacent 
sites creates an effective magnetic field. In the tight-binding approximation, the atomic 
motion is described by the Harper Hamiltonian, with a flux across each lattice plaquette 
that can realistically take any value between and n. We show how to take advantage of 
the superlattice to ensure that each plaquette acquires the same phase, thus simulating 
, a uniform magnetic field. We discuss the observable consequences of the artificial gauge 

field on non-interacting bosonic and fermionic gases. We also outline how the scheme 
£T) • can be generalized to non-Abelian gauge fields. 

> 

o 

1. Introduction 

The fractional quantum Hall (FQH) phases realized by two-dimensional electron gases 
in very large magnetic fields are among the most intriguing states of matter (see, for 
instance, [1]). In such systems, electrons "bind" to magnetic vortices, forming strongly 
correlated phases with striking properties, such as exotic excitations ("anyons") which 
obey fractional statistics (2]. Analogous quantum Hall phases should also arise in cold 
atomic gases when they are set into fast rotation (see [3] and references therein). Due to 
the mathematical similarity between Coriolis and Lorentz forces, rotating neutral gases 
are indeed the exact analogue of an assembly of charged particles plunged in a magnetic 
field. Observing these highly correlated phases is one of the major goals in the field of 
trapped quantum gases j3j H]. This goal has, however, not yet been reached, because 
of the difficulty of communicating the required amount of angular momentum (on the 
order of Nh per atom, with A" the number of particles) to the system j5j|6]. In practice, 
the residual static trap anisotropy limits the total angular momentum to much smaller 
values, for which the rotating gas is well described by a mean field approach [7] . 
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Recently alternative schemes to "simulate" artificial gauge fields for neutral atoms 
have been explored using two-dimensional (2D) optical lattices [H El [10J [TfJ [12], [T3| H%]. 
As they do not involve any mechanical rotation of the system, they should be less 
sensitive to the imperfection of the trapping potential and thus easier to implement. 
The guideline for these proposals is the celebrated Harper model [El [16], defined by the 
two-dimensional (2D) single-particle Hamiltonian, 



This tight-binding model was initially introduced [15] to describe electrons hopping on 
a square lattice perpendicular to a constant magnetic field B. The operator c\ m creates 
a particle at position (x = nd,y = md), d is the lattice spacing, and J is the tunnel 
energy to nearest neighbors in the absence of magnetic field. The effect of a vector 
potential (here in the Landau gauge, A = —Bye x ) is included in the so-called Peierls 
phase factor, 2nam, where a = eBd 2 /h is the flux per unit cell expressed in units of the 
flux quantum <p — h/e. The Harper Hamiltonian in Eq. (flj) has been studied extensively 
in the literature (see, e.g. [151 LTS1 [TF1 [IE])- For zero flux, the Harper model reduces to a 
standard tight-binding model on a square lattice, with a single Bloch band of width 8 J. 
For rational fluxes a = p/q (where p and q are integers), the model is still periodic but 
with larger unit cells of size qd x d. The Bloch band splits into q magnetic sub-bands, 
and the ground state becomes q— fold degenerate. This results in the peculiar self-similar 
structure in an energy- magnetic flux diagram known as Hofstadter's butterfly |16j . This 
structure exists in a strict sense for infinitely extended systems. For a finite system of 
size L, the details of this structure are washed out on fine scales such that q > L/d. 

In their seminal proposition for an implementation of ([T]) [8J, Jaksch and Zoller 
(JZ) suggested that the phase 2na can be imprinted by a laser beam inducing hopping 
between adjacent sites (see also [H [TUl HH [TJl [131 tH] f° r alternative proposals). A major 
advantage of this approach is that large flux a ~ 1 can be reached under experimentally 
realistic conditions. For a filling factor with ~ 1 atom per site, the quantum Hall regime 
then becomes reachable in presence of moderate atomic interactions [HI [191 I2Q1 [2lT [22] . 
Theoretical studies have found that some of these phases were essentially similar to that 
encountered in the bulk [HI [191 1201 [22] , whereas others are without counterpart in the 
bulk phase [231 El]. 

Although the existing proposals for implementing the Harper Hamiltonian with 
atomic quantum gases open many avenues, they are still challenging from an 
experimental point of view (for reasons discussed in detail below). In this paper, we 
extend the JZ proposal and discuss a new scheme based on an optical superlattice 
[251 [26] to generate a gauge potential leading to (JTJ. At variance with earlier works 
where alkali atoms were considered, we propose here to use atoms with a long-lived 
metastable excited state, such as alkaline-earth [271 [281 [29] or Ytterbium [301 [31] atoms. 
For concreteness we will discuss the case of Ytterbium, for which degenerate gases 
have already been produced for both bosonic and fermionic isotopes [301 EE]- We take 
advantage of this level structure to alleviate many practical difficulties encountered 
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Figure 1. (a) Sketch of the level scheme showing the two internal states g and e 
coherently coupled by a laser beam; (b): Dynamic polarizability a(X) of atoms in 
the ground (solid line) or excited states (dashed line). The light-shift experienced by 
the atoms is a(A)|E| 2 /2, with |E| the electric field strength at wavelength A. Two 
particular wavelengths are of special interest, a "magic" wavelength A m « 760 nm and 
an "anti-magic" wavelength A am rj 1.12 fim. (c) Sketch of the atomic configuration 
in a state-dependent optical lattice, showing the two interlaced sublattices for g (grey 
dots) and e (black dots) atoms. The y lattice is formed by a standing wave at the 
"magic" wavelength and confines both states identically. The x lattice is tuned near 
the "anti-magic" wavelength and confines the two states in two distinct sublattices 
shifted by A am /4. 



with alkalis. Importantly, our scheme uses only building blocks which have already 
been individually demonstrated. 

This paper is organized as follows. We first introduce the system, consisting of 
atoms with a ground state g and a long-lived metastable excited state e trapped in a two- 
dimensional (2D) lattice. The optical lattice is tuned near an "anti-magic" wavelength 
where the polarisabilities of e and g are opposite [32] . The two internal states thus live in 
spatially separated sublattices. We show that the case a = 1/2, leading to Dirac points 
analogous to those appearing in graphene [33], can be achieved rather simply with a 
single laser beam coupling e and g. We then extend our scheme to generate arbitrary 
values of a using an optical superlattice, at the expense of several lasers operating at 
different frequencies. Subsequently, we address the influence of the additional potential 
that must be added to confine the particles. We investigate its effect by looking at (i) 
the density of states that determines the density profile of a trapped ideal Fermi gas, 
(ii) the ground state of the system, relevant for non-interacting bosons. In both cases, 
we obtain a clear signature for a non-zero value of a. We conclude by indicating how 
the method can be generalized to generate non-Abelian gauge fields [34|. 

2. Two-electron atoms in state-dependent optical lattices 

Two-electron atoms, such as Ytterbium and alkaline-earth atoms, generally have a spin- 
singlet ( 1 S'o = g) ground state and an extremely long-lived spin triplet ( 3 Po = e) excited 
state (see Fig. The very long lifetime of e (around 20 s for Yb [35]) allows one to 
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operate optical atomic clocks at the resonance wavelength Ao ~ 578 nm between g and 
e [36] [37] 138] . It has also inspired novel proposals for quantum information processing 
[39| EOj HI] or quantum simulation [321 03]. States g and e usually react differently 
to far-off resonance laser light, as they couple optically to many different states. The 
polarizabilities of the two states calculated using the data from jS] are shown in Fig. [lb. 
For laser light tuned at the so-called magic wavelength A m (~ 760 nm for Yb), both 
states have the same polarizability and feel the same optical trapping potential [45J. 
Central to our proposal is the existence of an "anti-magic" wavelength A am ~ 1.12 /xm 
[32] . for which the polarizabilities of g and e are opposite J. As we will see, this allows 
to create easily internal-state dependent potentials. 

The first building block of our proposal consists in a laser beam at the magic 
wavelength, forming an horizontal light sheet that freezes the vertical motion of the 
atoms irrespective of their internal state e or g. From now on, we shall focus on the 
atom dynamics in the horizontal xy plane. The second ingredient consists in an optical 
standing wave along the y axis, also at the magic wavelength, which generates the same 
potential V(y) = —V cos 2 (ny / d y ) for both states g and e (V > 0). For the third 
building block, we apply another optical standing wave at the "anti-magic" wavelength 
A am along the x axis, generating the potential V'(x) = ±Vi cos 2 (Trx/2d x ), where the + 
and — signs refer to e and g, respectively. The resulting potential V'(x) + V(y) creates 
two shifted sublattices (see Fig. [lb). Atoms in state g sit at positions r g = (2n,m) and 
atoms in state e at positions r e = (2n + 1, m), where n, m are integers. Here we take 
as lengths scales d x = A am /4 along x, and d y = A m /2 along y. We write the Wannier 
functions for the lowest Bloch band as w g (r — r g ) and u> e (r — r e ). The potential V is 
chosen such that the tunnelling energy jw along the vertical lines is significant. On 
the opposite, we choose V± large enough to neglect direct tunnelling along the x axis. 

Tunnelling along the x axis is induced by a coherent optical coupling between g and 
e. Since atoms in e and g are trapped in separate optical lattices, this coupling induces 
hopping from one sublattice to the other. Following Jaksch and Zoller [8], we calculate 
the effective g — e hopping matrix element between two neighboring sites located at 
r g = (2n, m) and r e = r g + b = (2n + 1, m) as 

J(x) e iq.r 9 = ^iqT s J ^ _ h y^ Wg ^ ( 2 ) 

with Q and q the laser Rabi frequency and wavevector. We suppose that the coupling 
laser propagates in the yz plane, so that q • r g = In am, where |a| can take any value 
between and A m /2Ao (— 0.66 for Yb) by adjusting the angle between the coupling 
laser beam and the z axis. 

| There exists another anti-magic wavelength near 620 nm, which is unfortunately rather close to a 
transition connecting to the 3 Pq state (649 nm). A lattice at this wavelength is thus likely to suffer 
from excessive spontaneous emission. 

§ A many-electrons calculation 46 using more accurate parameters for the positions and widths of the 
optical transitions predict the anti-magic wavelengths near 1120.3 nm and 618.7 nm [A. Derevianko, 
private communication (2009)]. 
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Before proceeding with further analysis, let us give some numerical estimates for the 
different physical quantities involved in this experimental configuration. The anti-magic 
wavelength A am corresponds for Yb to a recoil energy E-^/h = h/(2m\ 2 m ) « 900 Hz. 
Realistic values up to several tens of kHz can be achieved for the depth of the optical 
lattice potential Vi/h. In order to avoid excitations to higher- lying Bloch bands, 
the Rabi frequency Q must be small compared to A/h, where A is the energy gap 
between the ground and first excited band. To give numerical estimates, we take 
V± = 20 En, corresponding to a gap A 8 Er. Assuming a moderate Rabi frequency 
Q = E-&/H, we find a laser-induced tunneling energy Jeg pa 0.05 E-&, much larger 
than the residual tunneling energy in each sublattice without change in internal state, 
Jgg = Jee^ ~ 2.5 x 10 _3 i?R. In principle, besides resonant transitions between the 
fundamental bands of the two sublattices, the coupling laser also induces second-order 
tunneling processes, where an atom tunnels to a neighboring site in a given sublattice, 
via a virtual transition to an intermediate excited band in the other sublattice. We 

(x) 

have verified that the corresponding tunneling rates were negligible compared to J gg 
for the range of parameters investigated here. Finally, one can calculate the spontaneous 
emission and associated heating rates for atoms in each sublattice [47J. For lattice depths 
V = Vi ~ 30 Er, we find total heating rates 7^ ~ 0.4 nK/s and 7^ ~ 0.9 nK/s, 
respectively, or 7^ / E R ~ 0.01 s _1 and 7^ /.Er ~ 0.02 s _1 . The heating rates 
(dominated by the vertical arm at A m ) are small enough to maintain the atoms in 
the ground band for several seconds. 



3. State-dependent lattice: Harper Hamiltonian for a — ~ 

We first consider the simplest experimental arrangement, with a spin-dependent lattice 
as described above and a single coupling laser (Fig. EK)- The two-component g — e 
quantum gas, which is assumed to populate the lowest Bloch band only, is described by 
the Hamiltonian 

it _ j \ y ( i2iramrf p, , i2vram-t A _i_ 1, n \ 

12 — J / j \ c c 2n+l,m c 2n,m ~r e, c 2n-l,m c 2n,m T U.l>. I 

n,m,n even 

+ h.c), (3) 

n,m,± 

when the lattice amplitudes are adjusted so that = j|g = J. As in [8], the 
operator & nm creates an atom at site (n, m), the parity of n identifying the internal 
state unambiguously. Alternatively, the internal degrees of freedom can be eliminated 
by going to a dressed state basis, and the spatial wavefunctions associated with each 
dressed state obey the same Harper equation [15l[T6] as the one derived from (pQ). Hence, 
at least in the absence of interactions (which we assume through this paper), the internal 
state is merely a label that can be ignored in the analysis. 

In the original equation (CTJ, the phase picked up across each link has the same 
sign for a given link direction in real space: e +t2narn in the +x direction and e~ i27Tam 
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in the — x direction. Here, we find in contrast that the phase is tied to the direction 
in internal space: e +t27Tam for g — > e transitions and e~ t27TOtm for e — > g transitions. In 
real space, given the sublattices geometry, this corresponds to a phase that alternates 
in sign from an arbitrary unit cell (or "plaquette" ) to the neighboring one along the x 
axis. More precisely, consider a particle going clockwise around an elementary cell (figure 
[lb). The phase picked up along the path (2n, m — l) g — > (2n,m) g — > (2n + l,m) e — > 
(2n+l,m-l) e (2n,m-l) g is the product of the two factors e ^^ e -^{m-\) =e i2na_ 
Thanks to the phase factor oc e jq ' rs , a particle making a loop around a plaquette acquires 
a non-zero global phase factor, reflecting the non-zero "magnetic flux" through the 
cell. However, for the adjacent plaquette, the path (also with clockwise orientation) 
(2n + 1, m - l) e ->■ (2n + 1, m) e -)> (2n + 2, m) s ->• (2n + 2, m - l) g -)■ (2n + 1, m - l) e 
leads to a phase factor e~* 27ra opposite to the previous one. One then achieves in this 
way a staggered magnetic field, with a phase per plaquette ±27ra constant on a given 
column and changing sign from one column to the next (Fig. [2b). 

The case a = | is peculiar, since in this case a staggered or a uniform flux both 
correspond to a tunneling phase e ±l2narn = (— l) m alternating between adjacent links 
depending on the parity of m. The Hamiltonian remains periodic, with a unit cell 
d x x 2d y and two sites per unit cell. The energy spectrum is given by 



where the Bloch vector (k x ,k y ) sits in the first Brillouin zone ] — Tr/d x , 7c/d x }x] — 
ir/2d y ,ir/2dy\. We note that this dispersion relation gives rise to two 'Dirac points" for 
(k x ,k y ) = (±7i/2d x ,Tc/2d y ) around which the dispersion relation is linear [T3], HH] (see 
also [IH [50j EU [52j [53] for discussions of the Dirac points occurring in different settings, 
such as hexagonal lattices or non-Abelian gauge fields). The physical consequences of 
the existence of these Dirac points are similar to those found in the case of Graphene 
[33J. However, their topological robustness [54] might be affected by small deviations 
from a = 1/2 that break the two-site periodicity. 

For a |, the staggered value of the phase per plaquette leads to a Hamiltonian 
that is different from ([T]), and which, in fact, does not exhibit any structure reminiscent 
of Hofstader's butterfly [55]. This can be linked to a general argument due to Haldane 
[56], emphasizing that breaking time- reversal invariance symmetry is the key factor to 
generate systems with properties similar to (integer) quantum Hall phases. The situation 
considered up to this point, where the same laser beam drives all g — > e transitions, is 
invariant with respect to time-reversal. We describe in the next section a mean to break 
this invariance and achieve a flux with a constant sign on all plaquettes. 

4. State-dependent superlattice: Harper Hamiltonian for arbitrary a 

4-1. Flux rectification 




(4) 



To ensure an identical phase 2ira per plaquette, instead of the staggered values ±2ira, 
we consider a situation where the on-site energies are modulated spatially along x with 
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Figure 2. Laser schemes for realizing staggered (a,b) and uniform (c,d) magnetic 
fields, (a): The figure shows the different on-site energies along x, in the case of 
a simple-period lattice. All transition frequencies are degenerate, so that a single 
laser beam at frequency ojl is sufficient to couple all sites across the lattice, (b): 
Staggered effective magnetic field for a single coupling laser propagating along y with 
degenerate transition frequencies as shown in (a), (c): In a configuration where an 
additional, doubly-periodic shift of the on-site energies is introduced, the degeneracy 
between the transition frequencies is lifted, (d): With several dedicated laser beams 
propagating along y, and uj'_ propagating along —y), one obtains a "rectified" 
uniform magnetic field. 



a period equal to twice the lattice spacing, as depicted in Fig. |2t (the y lattice remains 
identical for both states and is ignored in the following). Assume that the on-site 
energies are modulated according to 

E g — Vi, for x/d x = 4n, 
E e + V 2 , for x/d x = 4n + 1, 
Eg-Vx- V 2 , for x/d x =4n + 2, 
E e , for x/d x = 4n + 3. 

Here the E a (a = e, g) denote the internal energies in free space, where V\ is the 
amplitude of the fundamental lattice with period d x and where V 2 <C V\ is the amplitude 
of the modulation with period 2d x . In this potential landscape, the resonance frequencies 
for transitions linking neighboring sites become non-degenerate (Fig. Eb), 

lj + = cuq + (Vi + V-i)jTi for |<7;4n, m) — > |e; An + 1, m), 
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w_ = ujq + (Vi + 2V 2 )/h for |e;4n + l,m) -» \g; 4n + 2,m) , 
u + = ujq + (Vi + V<i)jh for |gr; An + 2, m) — > |e; 4n + 3, m), 
= w + Vi/h for |e; 4n + 3, m) — > \g; An + 4, m), 

with ItjjJq = E e — E g . 

To the state-dependent lattice, one thus applies three coupling lasers propagating 
along y with frequencies u + , U- and u'_. The laser at frequency u + is chosen 
with a wavevector k^, and the lasers at frequencies cu_ and u'_ with the opposite 
wavevector — k^. The alternance of the wavevectors compensates the alternance of the 
plaquette phases, thereby "rectifying" the staggered flux obtained in a single-frequency 
configuration. If we neglect off-resonant transitions, an analysis similar to that made in 
the previous section confirms that each plaquette now acquires the same phase factor 
e t2na , and that the system is indeed described by the Harper Hamiltonian (JT|). 

4-2. Optical superlattice 

To realize such a modulation in practice, we propose to superimpose along the x direction 
an additional lattice potential with a double spatial period 2A am , which can be written 
as V"{x) = Vg/ e cos(7rx/ 4d x +(p) 2 , where (p is the relative phase between the fundamental 
and double period lattices. Such superlattice potentials have been demonstrated 
experimentally for quantum gases [251 26J. This potential can be realized in practice by 
adding a second laser with frequency c/2A am phase locked to the fundamental lattice to 
control the relative phase <p [26]. The amplitudes V g / e depend on the internal state and 
are fixed by the initial choice of wavelength for the fundamental lattice. The analysis 
in Section 14.11 above (where V g = V e ) can be generalized straightforwardly to this less 
symmetric potential, with four non-degenerate transition frequencies in the general case. 
Moreover, it is possible to reduce the number of required laser frequencies to three by 
using a particular value of the relative phase <p that reproduces exactly the scheme 
described previously, with V 2 = V e cos ip — V g sin yjjj]. 

4-3. Conditions of validity 

We now discuss the domain of validity of our model. 

• First, as in the previous sections, we demand that laser- assisted tunneling dominates 
over regular tunneling without change in internal state. For a superlattice potential 
in the tight-binding regime, regular tunneling is suppressed when Vi > £7r, since 

Jgg i Jee <C -Er- 

• Second, the ratio MI/V2, with Q the Rabi frequency characterizing the atom-laser 
coupling, must remain small to neglect off-resonant transitions and describe the 

|| Practically, the fluctuations Sip of ip must be small enough to ensure that fluctuations of V2 remain 
much smaller than the width of the g — e resonance, set by the Rabi frequency ft, or Sip <C MI/V2 ^0.1. 
This is well within present technological capabilities. 
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system by an effective Harper model [8]. In practice, one can use a conservative 
value Ml = V 2 /10. 

• Finally, the modulation V2 must be small enough to avoid excitations to higher-lying 
Bloch bands by one of the coupling lasers. A reasonable choice is V2 = A/3, where 
A is the energy gap between the ground and first excited bands for V2 = 0. This 
guarantees that transitions mixing different bands are detuned by at least A/3. For 
example, if we focus on the |e;4n + 3,m) — > \g;^n + 4, m) transition, transitions 
linking the ground and first excited bands occur at frequencies ui'_±A/h~. One can 
check that the detunings for the different coupling lasers are | V2 ± A|, \2V2 ± A|, A, 
which are indeed > A/3 for the choice V% = A/3. 

We note that in the tight-binding regime, A is notably larger than which allows 
one to fulfill simultaneously all these conditions while maintaining a substantial value 
for J$. We summarize our considerations in Table [U where we give typical values 
for lattice parameters and laser-assisted tunneling matrix elements in the superlattice 
configuration. We note that the values for remain sizeable, being on the same 
order as the tunnel energies j'fg in the g sublattice (HQ = V2 = 0) at moderate depths 
(~ 10 E R ). 



Depth 


Energy gap 


Laser-assisted 




A/£ R 


tunneling Jig / £r 


10 


4.7 


16 x 10~ 3 


20 


7.8 


13 x HT 3 


30 


9.8 


9 x lO" 3 



Table 1. Characteristic tunnel energies (along the x axis) achievable for realistic 
experimental parameters in the superlattice configuration. Here J eg denotes the 
laser-assisted coupling from one sublattice to the other. We used a Rabi frequency 
Ml = V2/IQ = A/30 (see text). 

4-4- Comparison with the Jaksch-Zoller proposal 

The scheme as described in this paper is inspired by the earlier proposal by Jaksch 
and Zoller [8]. That proposal was designed for alkali atoms in spin-dependent optical 
lattices. In this context, g and e denoted two different hyperfine states in the ground 
state manifold, and the spin-dependent optical potential was obtained by exploiting 
the vector light-shift arising in a laser field with suitable polarization (see |57[ [58]. for 
example) . 

The first difference between our method and the JZ proposal comes from the choice 
of two-electron atoms, versus "one-electron" alkalis. To obtain a useful spin- dependent 
lattice for akalis demands in practice to tune the lattice laser between the D\ and D2 
lines, relatively close to resonance. For light alkalis, such as 6 Li or 40 K, the fine splitting 
between the D\ and D2 transitions is rather small (less than 1 nm), and spontaneous 
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emission rates are too large to maintain the cloud in the nK regime. The most practical 
candidate seems to be 87 Rb, which has only stable bosonic isotopes and would probably 
still suffer from spurious spontaneous emission due to the coupling (Raman) lasers or 
from the lattice lasers. The choice of Yb removes all these drawbacks, by using only far- 
off resonant lasers and a dissipation-free coupling transition, and by allowing to study 
bosonic and fermionic systems in the same setup. 

The use of a superlattice to ensure the same flux for all plaquettes constitutes 
the second main difference of our proposal with respect to JZ's proposal. There, a 
strong "tilting" potential independent of the internal state was introduced in order to 
distinguish \g; 2n, m) — > |e; 2n ± 1, m). Although this could be done in principle using a 
large static electric field or the light shift from an intense far-off resonant laser beam, it 
turns out to be rather challenging from an experimental point of vie In addition, this 
tilting potential has to be almost perfectly linear over the extent of the cloud. Otherwise, 
the various transition frequencies would depend on space, spoiling the flux uniformity 
and the scalability of the method. We believe that the use of a superlattice optical 
potential will be a more flexible and more reliable tool that can be readily implemented 
experimentally. 

5. Role of an external trap potential 

We now discuss the properties of a system of non-interacting particles described by 
Eq. ([I]) in order to identify observable consequences of the laser-induced gauge potential. 
For a uniform system without trapping potential, and for a = 0, the model reduces to a 
standard tight-binding Hamiltonian. The spectrum consists in an allowed energy band 
of width 8 J centered at zero energy. For the uniform Harper model, without trapping 
potential but with a = p/q ^ (p and q are integers), the initial Bloch band splits into 
q magnetic sub-bands separated by well-defined energy gaps. This results in q distinct 
peaks in the density of state (DoS-see Fig. (3^ and b). In a practical experiment, an 
additional trapping potential is inevitably present in addition to the optical lattice. 
Hence the physical system is both finite and inhomogeneous. This can strongly affect 
the properties of the uniform Harper model given in Eq. (H but also provide novel 
experimental signatures, as we discuss in this section. 

5.1. Density of states in presence of an external trap 

We consider here the Harper model in Eq. (CQ) with an additional isotropic trapping 
potential Vt rap = n(r/d) 2 . Let us recall that, with a = 0, one finds a quasi- continuum 

% The Stark shift from the electric field is proportional to a|i?| 2 , with a the static polarizability. Using 
a configuration where the cloud is offset from the center of the electrodes creating the field, one can 
create a linear potential of the form i]x/d x , with rj ~ aU 2 d x /D with D the distance between the 
electrodes and with U the applied electrostatic potential difference. For 87 Rb and D ~ 1 cm, we find 
that a site-to-site energy offset on the order of the energy gap A, ?y^A/3^10 kHz requires U ~ 50 kV, 
which is rather challenging to produce. 
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Energy [J] 



Figure 3. (a): Energy spectrum and (b): density of states (DoS) for the uniform 
model and a = 1/3. (c): density of states for the trapped model with a = 1/3 and 
confinement strength k/J = 0.01. The dashed line shows the expected density of 
states for E 3> J, irrespective of a. (d) Density profile for non-interacting fermions in 
a harmonic trap with k/J — 0.01, for a flux a = 1/3 (solid lines) and Fermi energies 
E-p = (half-filled band) and 4 J (filled band). The density profiles for a = (dashed 
lines) are also shown for comparison. The calculation was performed on a 63 x 63 space 
grid. 

of allowed states in the combined lattice plus external trap, in striking contrast with 
the uniform case [59]. States with energy lower than 4 J are qualitatively similar to the 
Bloch eigenstates of the uniform problem. These states are approximately contained 
in a disk of radius R e g ~ a/ J/k. States with higher energy (which corresponds to the 
bandgap of the uniform system) are well-localized states offset from the cloud center 
[60] . with an energy mostly determined by the local potential energy near the turning 
point. In two dimensions, the density of states saturates at an asymptotic value of tt/k. 

We now address the case with a = p/q ^ 0. The finiteness of the central region 
where tunneling is significant limits the relevant values of q to q < q max = R e s/d. We 
consider here the case a = 1/3 for concreteness. Solving the full model numerically, we 
find that the trapping potential has a dominant effect on the DoS and little remains 
of the initial peak structure (compare Fig. [3b and Fig. [3b). This indicates that global 
observables (such as thermodynamic quantities) will not be very sensitive to the presence 
of the gauge potential. We can understand this behavior from a semiclassical argument. 
Indeed, in a local density approximation, the density of states in a generic power-law 
potential n(r /d) 1 can be approximated as 
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where e n (k) = 4J"x n (k) is the relation dispersion for the nth magnetic sub-band of the 
untrapped system. This can be simplified into 

Here Sbz = {2it/d) 2 is the surface of the Brillouin zone. For a harmonic potential 
(7 = 2), the dependance on the reduced dispersion function \ is mild since it only arises 
from the upper bound of the integral: the peculiarities of the spectrum are essentially 
washed out by the spatial integration. This can be improved by using a steeper trap 
with 7 > 2, e.g. a quartic potential as demonstrated in [6]. 

5.2. Spatial distribution for non-interacting fermions 

In the case of non-interacting fermions in a gauge potential, a more sensitive observable 
is given by the particle density [22], assuming it can be measured in situ. Indeed, still 
within a semiclassical approximation, the spatial density can be written as 

n(r) = E / |^ /fd ^ + WO - Ef] , 

with /fd the Fermi-Dirac distribution at T = and E F the Fermi energy. We rewrite 
this as 

d (2) k ,00 



n r 



— jjE / FD [V trap + E - E F ) 6 [e n (k) - E) 

dE'f FD [Krap + ^' - ^ F ] Phom [E'\ , (7) 



n 
oo 



1.7 



where 



^ = Ey j2^n^k)- E ] (8) 



Pho 



is the density of state for the homogeneous system. For low enough temperatures, the 
clustering of energy levels in subbands in presence of a gauge potential will result in well- 
defined plateaux in the spatial density of the trapped cloud. These plateaux correspond 
to "band-insulating" regions, similar to that appearing for other incompressible systems 
in traps such as Mott insulators J6TJ [62]. We have verified this fact in our numerical 
calculations, as shown in Fig. (3U. We conclude that the in-trap density profile will 
provide a clear signature of the gauge potential in the case of fermionic atoms. 



5.3. Momentum distribution for non-interacting bosons 

With non- interacting bosons, such a signature will be absent, as the atoms will condense 
into the lowest single-particle states (for the purpose of this paper, we do not consider 
interactions). However, signatures of the gauge potential will appear in reciprocal 
momentum space, which can be probed using standard time-of-flight imaging. The 
dephasing term e l2lTam generates quasi-momentum components at the various harmonics 
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Figure 4. Typical momentum (a) and spatial (b) distributions in the ground state 
for non-interacting bosons and a flux a — 1/3. The presence of a gauge potential 
reveals itself in momentum components at ±27r/3<i in the former plot, and in the fine 
structures with period 3c? in the latter. The confinement strength is k/J — 0.005. 

of the fundamental period 2ira. As a result, clear satellite peaks can be observed in the 
momentum distribution (see Fig. Hk), which will provide a direct experimental signature 
for bosonic species. We note in addition that the density distributions corresponding 
to these eigenstates shows periodic structures with period q (see Fig. 0b), although this 
signature appears on a much finer scale than in the fermionic case (thus being more 
difficult to observe). 

6. Conclusion and perspectives 

In this article, we have discussed how to implement Abelian artificial gauge potentials for 
Yb (or alkaline-earth) atoms, a scheme valid for bosons and fermions alike. In contrast 
to the bosonic isotopes with zero total spin, fermionic isotopes have nuclear spin. This 
potentially enables to simulate non-Abelian gauge fields, along the lines proposed by 
Osterloh et al. [34). Let us focus for concreteness on the 171 Yb isotope, with nuclear 
spin I = |. In the presence of a moderate magnetic field of a few tens of Gauss, the 
degeneracy of the mj = ±| states within the e and g manifolds are lifted, giving in 
general four different internal transitions, ±~ — > ±~ and ±~ — > =f|, with different 
frequencies. At the expense of using more laser frequencies and state-dependent lattices 
along both x and y directions, one can then arrange for the lasers to imprint different 
phases depending on the internal state or generate arbitrary rotations in internal space. 
As shown in [33], a straighforward generalization of the Abelian scheme then allows 
to engineer non-Abelian gauge potentials. Instead of a simple phase, hopping from 
sublattice g to e then corresponds to a rotation in internal space generated by non- 
commuting matrices. Combined with interatomic interactions, an entirely new class 
of superfluid or strongly correlated systems becomes accessible with ultracold atoms 
[331 E3J E31 E31 E5]. In a wider context, this technique can be used to emulate spin-orbit 
coupling, similar to the Rashba interaction oc a x ■ k y — a y ■ k x in planar semiconductors, 
or for the study of strongly interacting, ultracold fermionic matter in a non-Abelian 
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background field, possibly relevant to high-energy physics. 
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